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ABSTRACT 


We present a class of compensators in explicit form (not requiring numerical computer 
calculations) for stabilizing flexible structures with collocated rate sensors. They are based on the 
explicit solution, valid for both Continuum and FEM Models, of the LQG problem for minimizing 
mean square rate. They are robust with respect to system stability (will not destabilize modes even 
with mismatch of parameters), can be instrumented in state space form suitable for digital 
controllers, and can be specified directly from the structure modes and mode “signature” 
(displacement vectors at sensor locations). Some simulation results are presented for the NASA 
LaRC Phase-Zero Evolutionary Model — a modal Truth model with 86 modes — showing 
damping ratios attainable as a function of compensator design parameters and complexity. 

1. INTRODUCTION 


In this paper we present a class of compensators for stabilizing flexible structures with 
collocated rate sensors for Continuum as well as Finite Element or truncated Modal models. They 
are derived by solving explicitly the optimal control corresponding to an LQG problem. The 
Compensator Transfer Functions are strictly positive real and as a result they are robust with 
respect to system stability. They can be determined based solely on the mode frequencies and 
mode “signatures” (displacements at the sensor sites) and can be instrumented directly in “state- 
space” form. 


We begin in Section 2 with the LQG problem and its solution. Section 3 highlights the 
features of the Compensator Transfer Function. Section 4 is devoted to Continuum models where 
the transfer function is nonrational. The main result on controller design is in Section 5 which 
shows how to design the compensator from a modal model of the structure, and in particular, how 
to construct a hierarchy of compensators of increasing order. Simulation results, confined to 
stability properties (damping ratios), are by no means exhaustive and are presented in Section 6 based 
on the modal model of the NASA LaRC CSI Phase-Zero Evolutionary Model [4]I Noise response 
performance is not included but is expected to be good because of the LQG criterion optimization. 
Conclusions are in Section 7. 


* Research supported in part under NAS- 1-191 58, NASA LaRC. 
I References 1-6 are cited in text. 


475 



Foi lecent related work on controller design for collocated sensors see [5, 6| and the 
references therein. 


2. THE LQG PROBLEM AND ITS SOLUTION 


To state the LQG problem, we begin with the canonical time-domain dynamics of a flexible 
structure with collocated rate sensors which, whether it is a Finite Element Model (and hence finite 
dimensional) or Continuum Model (and hence infinite-dimensional) can be expressed in the form: 


Mx(t) + Ax{t) + Bu{t) + BN a (t) = 0 
v(f) = B*x(t) + N r (t) 


(2.1) 


where in the case of FEM, 

M is the mass matrix (nonsingular, nonnegative definite) 

A is the stiffness matrix (nonsingular, nonnegative definite) 

B is the control matrix (rectangular matrix) 

u(-) is the control vector (nx 1, assuming n actuators) 

■*(') is the “displacement” vector 

N a (-) is the actuator noise assumed white Gaussian with spectral 
density d a I , / being the nxn Identity matrix 

v(-) is the sensor output 

B* represents the transpose of B 

N r (-) is the sensor noise assumed white Gaussian with spectral density d r l . 


For the Continuum Model such a representation is still possible with *(•) now allowed to 
range in a Hilbert space (however complicated the structure), with A, M, and B being linear operators: 

M bounded linear, self-adjoint, nonnegative definite with M ~ 1 bounded; 

A closed linear, self-adjoint, nonnegative definite with compact resolvent, 

resolvent, the resolvent set including zero 

B maps E n Euclidean n-space into H, and 

B* represents the adjoint. 

See m, [2]. 


The LQG problem we shall consider is that of finding the control /<(•) (or compensator) that 
minimizes the mean square time average: 


476 



lim -Uj \\B*x(t)\\ 2 dt + MOW 2 dt\ ( 2 

7’-»00 l 1 Q Q J 

where ^ > 0. The optimal compensator transfer function (nxn matrix function) can be given in 
explicit form (see [1]): 

\| f(p) = gpB*(p 2 M + A + ypBB*) B , Re. p > 0 (2 


where 


U a / d r 


<X ’ 


Y = 


under the assumption that 


B*<$> k * 0 


for any k, where § k are the modes orthonormalized with respect to the mass matrix: 


A§ k — > 


m k , <t>*i = 1 • 


In the finite-dimensional (e.g., FEM) case, the compensator can be realized in the (finite- 
dimensional) state space form: 

u(t) = gB*Y(t) 

MY(t) + AY(t) + yBB*Y(t) = Bv(t) . 

The corresponding mean square control power: 


l im U IHOII 2 A = ^ Jr. . 


The corresponding mean square displacement: 


i 7 

lim ^ J \\B*x(t)\\ 2 dt 


tijd, + ^Mtx.B^B. 


Formulas (2.9) and (2.10) hold as well in the infinite-dimensional (Continuum Model) case. 


3. FEATURES OF THE COMPENSATOR TRANSFER FUNCTION 


Some significant features of the compensator transfer function which are noteworthy are: 

(i) As ^->0, \|/(p) -)Nd a ld r we note that 

Mp) = f ^ d a ld r 

is the optimal “static” or “direct connection” or “PID” controller. Note that as X -> 0, the control 
power given by (2.9) becomes infinite, as we expect. 
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(ii) \j/(p) is “positive real” — that is to say: 

\|/(p) holomorphic in Re.p > 0 

V(/>) + V(P)* nonsingular, and positive definite, for Re. p > 0 
where * denotes conjugate transpose. Of course \j/(-) e K x . 

For the importance of positive realness for robustness, see [6]. 

Let a compensator transfer function be defined by (2.3) where g and y are arbitrary subject 
only to the condition that 

(3.1) 

Then for d a = 0 (no actuator noise), the corresponding mean square displacement 

T 

T / \\B*x(?)\?dt = ^Tr .(fiM-'fi). (3 . 2 ) 

This follows from [ 1, (6.11)]. 

4. EXAMPLES: CONTINUUM MODELS 


For a given Continuum Model for structures (see [1], [3]) it is possible to reduce (2.3) further 
to yield finite-dimensional matrix transfer functions which are, however, not rational functions of p 
Thus (2.3) becomes: 

V(p) = gpB*(p 2 M b +T(p) + ypB u B*)~'B u 

where 

T(p) . self-adjoint nonnegative definite nonsingular matrix for p > 0. 
Nonrational meromorphic function of p. Re p > 0. 

B u is d x n where 

d : number of nodes x dimension of displacement vector at each node 

n : number of actuators/sensors 

M b is the dxd mass/moment of inertia matrix corresponding to the nodes. For the SCOLE 
configuration [1], d = n = 5, 

R R* - J 
D u°u ~ *5x5 

and an explicit expression for T{p) is given in [3], and involves hyperbolic sine and cosine 
functions of p. An example where d* n is given by the NASA LaRC Phase Zero Evolutionary 
Model [41 for which n = 8 and d = 48. See also [3] for some “textbook” examples for d = n = 1 , 
Bernoulli beam bending or torsion. 

The system modes are given by 

M bP 2 + T(p) = 0 . 
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5. MODAL APPROXIMATION 


Let {<))* } denote the eigenvectors (mode “shapes”) of the stiffness operator A with co* the 
corresponding eigenvalue (angular frequency). Let b k denote the corresponding mode “signature” 
row vector: 

b k = (B*<t>*)* lxn row vector 


where n is the number of actuators. Let for each N > 1 , 


b i 


B n - 



D n = Diagonal (oof, ..., C0yv) . 


For arbitrary g,y>0, define the compensator transfer function 

Vn(p) = gP B tS{P 2 h + d n + 1 P b n b $Y' b n > Re. p > 0 . 


(5.1) 


Then 


(a) \|/ N (-) is positive real as soon as is nonsingular; 

(b) For g, y defined by (2.4), \\f N (p) converges to the optimal compensator \|t(p) 
given by (2.3) as N ->oo (and holds a fortiori in the finite-dimensional case, 
where the sequence terminates). 


Note that the modal approximation requires only the modes and the “modal signature”: 
modal displacement at the sensor sites. Note also that (5.1) automatically yields a strictly positive 
real rational transfer function approximation for the case of the Continuum Model — yielding, in 
fact, a new technique for such approximation. Moreover it has the direct state space representation: 


u(t) = gBfiY(t) 

Y(t) + D N Y(t) + y B N BftY(t) = B N v(t) . 


(5.2) 


It is also important to note that for a given finite-dimensional modal model with m modes, 
say, we can choose any N modes for the approximation, not necessarily the first N . Moreover the 
stability properties of the system are determined by the (“closed-loop”) eigenvalues of the “system” 
2 (m + N) x 2 (m + N ) matrix: 


xm 

I m Xm 

OmxN 

®mxN 

-Dm 

0/n Xm 

^m xN 

-g B m B $ 

®Nxm 

fy/xm 

®NxN 

InxN 

®Nxm 

b nK 

-D n 

~J b n b n 
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This is readily seen to be a stable matrix under our condition that B N Bft is nonsingular and 
g, y > 0. Moreover, including a damping matrix D > 0 in the Truth Model (replace 0 mXm in 
(5.3) by -D), we see that 

Trace W = -Tr.D - y Tr. BftB N = 2(sum of real parts of eigenvalues) (5.4) 

again illustrating the robustness. Also we have 

product of roots = \W\ = |D m | • |D W | (5.5) 

where |*| denotes determinant. Finally let us note that the eigenvalues are the roots of 

\pl - W\ = | p 2 +D m + pD\ | p 2 + D n + y pB N Bj$ + P 2 gB N B * ( p 2 + D m +pDy 1 B m Bft | = 0 (5.6) 

where the first factor is an mxm determinant and the second factor is an NxN determinant, if 
g = 0 structure modes are unaffected. We can see from (5.4) - (5.6) that the total damping increases 
as y is increased but the damping in the structure modes can decrease depending on how large y is. 
Thus, for each fixed value of the gain g, there is apparently an optimal choice for y, which may 
depend on the mode frequency in general. 

6. SIMULATION RESULTS: PHASE ZERO EVOLUTIONARY MODEL, NASA LaRC 


For evaluating control performance by simulation, the NASA LaRC CSI Phase-Zero 
Evolutionary model data [4] is used — specifically, the modes and modal signatures. 

We denote by m the number of modes in the truth model of the structure, and let N = the 
number of modes in the compensator, as in (5.2) and (5.3). The compensator is thus characterized by 
the “gain” parameter g, the “damping” parameter y, and N the number of “control” modes 
(or 2 N = number of states). Note that for each N 


Vn (p) kl ( ~ “direct connection” or “static controller”) 
as g -t<», y -too, keeping the ratio 

8 


(6.1) 


fixed. Also for N = m, we can use (3.2) to determine the mean square displacement where now 

m itij j 1 2 

Tr. B*A~ 1 B = I — j , 


so that the mean square displacement can be expressed 

d J± V II** H 2 

2y \ (4 ’ 

and thus increases with g and decreases as y increases. 


Y 2 > 4 g , 


( 6 . 2 ) 
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Data for the modal model of the NASA LaRC Phase Zero Evolutionary Model is taken from 
[4]. Here n = 8 so that each b k is 8x1 , and m for the Truth model is 86, and D = 0 (no 
damping). In this case the optimal compensator will use all 86 modes. Table I shows the mode 
frequencies and a typical mode shape (b-j). is nonsingular. 


Figures 1, 2, and 3 show the behavior of the damping ratio for fixed g as a function of the 
damping parameter y, for angular frequencies 314 (the 86th mode), 106 (the 43rd mode) and 9.25 
(the 7th mode), respectively. Note the occurrence of the maximum for all around y = 8. Figures 4, 
5, and 6 show the attainable damping ratios for a compensator with N = 30, and the same gam 


§ Personal communication, S. M. Joshi, NASA LaRC, 1992. 


Mode 

Number 

1 = 

2 - 

3 = 

4 = 

5 = 

6 = 

7 - 

8 = 

9 = 

10 = 
11 = 
12 = 

13 = 

14 = 

15 = 

16 = 

17 = 

18 = 

19 = 

20 = 
21 = 
22 - 


(Angular) 

Frequency 

Mode 

(Angular) 

Frequency 

rad/sec 

Number 

rad/sec 

0.9243 

23 = 

38.8269 

0.9365 

24 = 

39.1489 

0.9752 

25 = 

40.6580 

4.5872 

26 = 

41.9080 

4.6985 

27 = 

46.3219 

5.4913 

28 = 

52.1080 

9.2580 

29 = 

52.8380 

10.9209 

30 - 

53.1319 

11.8310 

31 = 

55.4410 

14.4600 

32 = 

56.0859 

15.9259 

33 = 

56.3370 

17.8349 

34 - 

58.0239 

21.4850 

35 = 

59.8600 

21.9050 

36 - 

62.2190 

22.5400 

37 = 

78.4509 

25.2250 

38 = 

85.5439 

25.3349 

39 = 

89.9440 

26.4249 

40 = 

92.4779 

27.5949 

41 = 

99.8170 

31.6009 

42 = 

105.8899 

31.6289 

43 = 

106.7699 

34.5660 

44 = 

116.1500 


Mode 

(Angular) 

Frequency 

Mode 

(Angular) 

Frequency 

Number 

rad/sec 

Number 

rad/sec 

45 = 

120.5899 

67 = 

195.1799 

46 = 

133.1399 

68 = 

196.9100 

47 = 

137.8399 

69 = 

198.1300 

48 = 

139.5700 

70 = 

198.7400 

49 = 

147.2599 

71 = 

203.3300 

50 = 

154.2500 

72 = 

209.7299 

51 = 

156.5399 

73 = 

231.7400 

52 = 

156.9199 

74 = 

232.9499 

53 = 

161.0299 

75 = 

239.8899 

54 = 

164.6499 

76 = 

241.4499 

55 = 

165.1399 

77 = 

244.9299 

56 - 

166.9100 

78 - 

244.9900 

57 = 

170.8300 

79 = 

247.8300 

58 = 

173.7700 

80 = 

255.7700 

59 = 

180.0099 

81 = 

268.6499 

60 = 

181.9100 

82 = 

270.5499 

61 = 

183.1300 

83 = 

275.7200 

62 = 

186.6199 

84 = 

287.5299 

63 = 

187.7500 

85 = 

308.7500 

64 = 

191.2400 

86 = 

314.2399 

65 = 

66 = 

192.0599 

193.7899 




Mode Shape b n 


- 0.9084100127220 

0.0009075700073 


- 0.0009229300194 

- 0.4345000088215 


0.4756200015545 

- 0.0070441002026 


0.0040255999193 

1.1181999444962 


Table I. Mode Frequencies And Sample Mode Shape 
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g ( 20), but here there is no maximum for the mode frequency of 314 which is not included in 

the controller modes. For the mode frequencies 9 and 106, the maximum occurs around y = 5 and 
y = 30, respectively. Table II shows the damping ratios for all (angular) frequencies (both control 
and structure modes; the former have higher damping) between 1 and 150 for g = 20 and y = 5 for 
the optimal compensator ( N = 86). Table III shows the same for N = 30. The damping ratios are 
seen to compare favorably with those reported in [5, 6], depending of course on the appropriate gain 
setting, but detailed comparative evaluation will need further study. Table IV shows the damping 
ratios for zero gain ( g = 0) with N = 86 and y = 10 which should help distinguish the control 
modes from structure modes (the latter have zero damping). Figure 7 shows the damping ratio as 
a function of g for N = 86, y = 10 and mode frequency 21.485. The damping increases for the 
structure mode and decreases for the control mode. 


10*5 .Damping Ratio 



Gamma 


Figure 1. Damping ratio vs. damping parameter y: mode angular frequency 314; N = 86. 
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Damping Ratio 



Figure 6. Damping ratio vs. damping parameter y: mode angular frequency 9.25; N = 30. 


10*4 x Damping Ratio 



Figure 7. Damping ratio vs. g: mode frequency 21.49 r/s; N = 86; y = 10; C: control; S: structure 
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Angular Damping 
Frequency Ratio 

. 1 7 1 06E+01 .4881 4E+G0 
. 2 52 33E+01 . 4234 9E+0Q 
. 38285E+01 . 18 903E+00 
. 40291E+01 .S3048E+00 
.50981E+01 . 53701E-0 1 
. 550 10E+01 . 27800E+00 
. 62846E+01 . 4 8 991E+00 
. 65389E+01 . 15363E+00 
. 735 98E + 01 . 25791E+00 
. 88473E+01 . 2857 6E+00 
. 96096E+01 . 23182E+00 
. 98633E+01 . 46814E-01 
. 1 4 333E+02 . 26440E-02 
. 1 4 370E+02 . 13551E-01 
. 1 5912E+02 . 94853E-03 
. 1 5918E+02 . 12789E-03 
. 1 68 58E+02 . 34 110E+00 
. 1 74 38E+02 . 154 97E-01 
. 17615E+02 . 364 68E-02 
. 18676E+02 . 2751 8E+00 
. 2 1 4 78E+02 . 2698 8E-03 
. 214 91E+02 . 25982E-03 
. 21897E+02 . 29663E-03 
. 2 191 3E+02 . 24 238E-03 
. 22536E+02 .20344E-03 
. 22544E+02 .20967E-03 
. 23354E+02 .22881E+00 
. 23666E+02 .31516E-01 
. 25234E+02 .10998E-02 
. 25241E+02 . 24 G57E-03 
. 25492E+02 .34385E-01 
. 2 6321E+02 . 16274E-02 
. 26361E+02 . 82509E-03 


Angular Damping 
Frequency Ratio 

, 27975E+02 . 10362E+00 
. 30735E+02 .22514E+00 
. 30903E+02 . 16990E+00 
. 31554E+02 . 10862E-02 
. 31573E+02 . 34084E-02 
. 31601E+Q2 . 36615E-05 
. 31601E+02 . 6047 5E-05 
. 332 65E+02 .29600E-01 
. 35339E+02 .12758E-01 
. 362 45E+02 .22248E-01 
.38097E+02 .15841E-01 
. 394 85E+02 .23587E-01 
. 397 15E+Q2 . 64 4 62E-02 
. 40253E+02 . 96291E-02 
. 4 2723E+02 . 40264E-01 
. 47595E+02 .20357E+00 
.491 16E+02 . 37251E-01 
. 50628E+02 .74785E-01 
. 52567E+02 .48805E-03 
. 52612E+02 . 36127E-03 
. 52923E+02 . 23978E-01 
. 53030E+02 . 65234E-03 
. 53116E+02 . 39229E-02 
. 55388E+02 .51163E-03 
. 554 76E+02 . 11573E-02 
. 55627E+02 .47474E-01 
. 55777E+02 . 20226E+00 
. 56012E+02 . 50958E-03 
. 56120E+02 . 32213E-02 
. 57973E+02 .33939E-03 
. 58058E+02 .22112E-02 
. 59857E+02 .16844E-04 
. 59860E+02 .10358E-03 


Angular 

Damping 

Frequency 

Ratio 

. 61 950E+Q2 

. 10323E + 00 

• 62123E + 02 

. 77622E-03 

. 62301E + 02 

. 15341E-02 

• 77219E + 02 

. 89300E-02 

. 79869E+02 

. 16028E-01 

. 85536E+02 

-45737E-04 

. 855S4E+02 

. 13440E-0 3 

- 89610E+02 

. 18167E-02 

. 90357E+02 

. 56712E-02 

. 92127E+02 

. 21810E-02 

. 92875E+02 

. 3577 5E-02 

. 99737E+02 

. 42099E-03 

. 99912E+02 

. 10766E-02 

. 10474E+03 

. 624 17E-02 

. 10629E+03 

. 21788E-02 

. 10703E+03 

. 53890E-02 

. 10757E+03 

. 10741E-01 

. 11586E+03 

. 14946E-02 

. 11648E+03 

.24414E-02 

. 12059E+03 

. 22062E-04 

. 12059E+03 

.32178E-04 

. 13306E+03 

. 38080E-03 

, 13322E+03 

. 42566E-03 

13759E+03 

. 12 103E-02 

13809E+03 

. 12406E-02 

13957E+03 

. 39718E-05 

13957E+03 

. 50153E-05 

14684E+03 

. 17338E-02 

14771E+03 

. 22478E-02 


Table II. Damping Ratio Vs. Angular Frequency: N = 86; g = 20; y = 5 


Angular 

Damping 

Angular 

Frequency 

Ratio 

Frequency 

. 17111E+01 

. 4 8 787E + 00 

. 27681E+02 

. 2 52 40E+0 1 

. 42289E + 00 

. 30664E+02 

. 38283E+01 

. 18850E + 00 

. 30908E+02 

40268E+01 

. 53155E+00 

. 31554E+02 

50982E+01 

. 53764E-01 

.31573E+02 

55074E+01 

. 2764 6E+00 

. 31601E+02 

62704E+01 

. 49202E+00 

. 31601E+02 

6541 3E+01 

. 15327E+00 

. 332 68E+02 

73966E+01 

. 2558 9E+00 

. 35300E+02 

88395E+01 

. 28765E+00 

. 3624 6E+02 

961 14E+01 

. 2420 4E+QQ 

. 380 90E+02 

98657E+01 

. 46669E-01 

. 39472E+02 

14333E+02 

.26 268E-02 

. 397 16E+02 

14371E+02 

. 13622E-01 

. 40247E+02 

1 59 12E+02 

. 94 628E-03 

. 42747E+-02 

1 59 18E+02 

. 12646E-03 

. 4 7582E+02 

168 7 3E+02 

. 34829E+00 

. 49299E+02 

17439E+02 

. 15483E-01 

. 50456E+02 

17615E+02 

. 36319E-02 

. 52567E+02 

18473E+02 

. 29242E+00 

. 52615E+02 

21478E+02 

. 26755E-03 

. 53031E+02 

21491E+02 

. 26123E-03 

. 53108E+02 

21897E+02 

. 29369E-03 

. 554 46E+02 

2 1913E+Q2 

. 24459E-03 

. 55694E+02 

22536E+02 

. 19617E-03 

. 55773E+02 

22544E+02 

.2124 5E-03 

. 561 12E+02 

23367E+02 

. 23050E+00 

. 58052E+02 

23680E+02 

. 29890E-01 

. 584 17E+02 

25234E+02 

. 11092E-02 

. 598 61E+02 

2524 1E+02 

. 24 12 5E-03 

. 62235E+02 

2 54 95E+02 

. 33 94 6E-01 

. 78798E+02 

26323E+0 2 

. 1 6797E-02 

. 85547E+02 

26362E+02 

. 80854E-03 

. 90088E+02 


Damping 

Angular 

Damping 

Ratio 

Frequency 

Ratio 

. 10975E+00 

. 92550E+02 

. 19914E-03 

, 22748E+00 

. 99846E+02 

. 13752E-03 

, 16861E+00 

. 10612E+03 

. 41269E-03 

10567E-02 

. 1068 9E+03 

. 2357 IE -03 

34060E-02 

. 11618E+03 

. 42708E-04 

36214E-05 

. 12059E+03 

. 81852E-06 

60906E-05 

. 13315E+03 

. 42 638E-05 

29376E-01 

. 13785E+03 

. 64145E-05 

13261E-01 

.13957E+03 

. 88 491E-07 

21919E-01 

. 14728E+03 

. 89389E-05 


.15162E-01 
. 24008E-01 
. 61748E-02 
. 98628E-02 
. 38763E-01 
. 20393E+00 
. 30838E-01 
. 78217E-01 
. 5074 6E-03 
. 32023E-03 
. 60754E-03 
. 36922E-02 
. 30132E-03 
. 20493E+00 
. 44503E-01 
. 13185E-02 
. 7257 3E-03 
. 24 337E-01 
. 78116E-04 
. 19700E-03 
. 19928E-02 
. 18896E-04 
. 74 096E-03 


Table III. Damping Ratio vs. Angular Frequency: N = 30; g = 20; y = 5 
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Angular 

Frequency 

. 45872E+01 
. 46986E+01 
. 53549E+01 
. 54 914E + 01 
. 89044E+01 
. 9258 1E+01 
. 10921E+02 
. 1 1831E+02 
. 13632E+02 
. 14263E+02 
. 14460E+02 
. 15914E+02 
. 15926E+02 
. 1 7552E+02 
. 17835E+02 
.214 85E+02 
.21487E+02 
. 21905E+02 
. 21907E+02 
. 22375E+02 
. 22540E+02 
. 22541E+02 
. 24 4 59E + 02 
. 25225E+02 
. 2 5243E+02 
. 25335E+02 
. 2 6335E+02 
. 2 6425E+02 
. 2 7595E+02 
. 304 56E+02 
. 315 63E+02 
. 31601E+02 
. 31601E+02 


Damping 

Ratio 

. 0000 0E-»-00 
. OOOOOE+OO 
. 10 629E+00 
. OOOOOE+OO 
. 8 4 2 4 9E-01 
. OOOOOE+OO 
. 00000E+00 
.OOOOOE+OO 
.17079E+01 
. 77 62 9E-02 
. OOOOOE+OO 
. 28005E-03 
. OOOOOE+OO 
. 82357E-02 
. OOOOOE+OO 
. OOOOOE+OO 
. 78458E-03 
.OOOOOE+OO 
. 90374E-03 
. 13007E+00 
.OOOOOE+OO 
. 24 691E-03 
. 13332E+00 
.OOOOOE+OO 
. 45734E-03 
. OOOOOE+OO 
. 13181E-02 
. OOOOOE+OO 
. OOOOOE+OO 
. 13725E+00 
. 3054 4E-02 
. OOOOOE+OO 
. 11152E-04 


Angular 
Frequency 
. 31629E+02 
. 34566E+02 
. 352 59E+02 
. 387 62E+02 
. 38827E+02 
. 39149E+02 
. 397 44E+02 
. 40658E+02 
. 4 11 18E+02 
. 4 19 08E + 02 
. 4 6322E + 02 
. 4 66 30E+02 
. 50739E+02 
. 52108E+02 
. 52588E+02 
. 528 38E+02 
. 52 9 88E+02 
. 53132E+02 
. 55387E+02 
.554 41E+02 
. 55964E+02 
. 56086E+02 
. 56337E+02 
.57938E+02 
. 58024E+02 
. 59855E+02 
. 598 60E+02 
. 621 25E+02 
. 622 11E+02 
. 77305E+02 
. 78451E+02 
. 85533E+02 
. 8554 4E+02 


Damping 

Ratio 

. OOOOOE+OO 
. OOOOOE+OO 
. 24 327E-01 
. 32907E-01 
.OOOOOE+OO 
. OOOOOE+OO 
. 22024E-01 
. OOOOOE+OO 
. 24 327E+00 
. OOOOOE+OO 
.OOOOOE+OO 
. 20336E+00 
. 80977E-01 
.OOOOOE+OO 
. 42582E-03 
.OOOOOE+OO 
. 21608E-02 
.OOOOOE+OO 
. 22200E-02 
.OOOOOE+OO 
. 17243E-02 
.OOOOOE+OO 
.OOOOOE+OO 
. 11515E-02 
.OOOOOE+OO 
. 51996E-Q4 
.OOOOOE+OO 
. 34738E-02 
.OOOOOE+OO 
. 39799E-0 1 
.OOOOOE+OO 
. 19002E-03 
. OOOOOE+OO 


Angular 
Frequency 
. 89479E+02 
. 89944E+02 
. 921 85E+02 
. 92478E+02 
. 99727E+02 
. 998 17E+02 
. 10466E+03 
. 10589E+03 
. 10616E+03 
. 1 0677E+03 
. 11595E+03 
. 1 1615E+03 
. 12059E+03 
. 12059E+03 
. 13312E+03 
. 13314E+03 
. 13783E+03 
. 137 84E+03 
. 13957E+03 
. 13957E+03 
. 14705E+03 
. 14726E+03 


Damping 

Ratio 

. 73662E-02 
.OOOOOE+OO 
. 99239E-02 
.OOOOOE+OO 
.17 4 4 3E-02 
.OOOOOE+OO 
. 29513E-01 
.OOOOOE+OO 
. 92339E-02 
.OOOOOE+OO 
. 63634E-02 
.OOOOOE+OO 
. 97457E-04 
. OOOOOE+OO 
. 15 99 7E-02 
. OOOOOE+OO 
. 49739E-02 
.OOOOOE+OO 
. 16752E-04 
.OOOOOE+OO 
. 77790E-02 
.OOOOOE+OO 


Table IV. Damping: Zero Gain: V=86; y= 10 


7. CONCLUSIONS 


A class of compensators in explicit form (not requiring computer calculations) has been 
presented for stabilizing flexibility structures with collocated rate sensors. They are optimized for 
the LQG criterion for minimizing the mean square rate and hence have inherently good noise 
response features. They are robust with respect to system stability, can be instrumented in state 
space form suitable for digital control and above all can be specified to any complexity desired 
directly from the structure modes and mode signatures at the sensor sites. Simulation results are 
presented for the modal model of the NASA LaRC Phase Zero Evolutionary Model — mainly 
damping ratios attainable and their dependence on compensator design parameters and complexity. 
The damping ratios compare favorably with those reported in [5, 6], but any detailed comparative 
evaluation of course is possible only after further study. 
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